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Abstract 



^ \ A generalization of the factorization technique is shown to be a powerful 

^^ ' algebraic tool to discover further properties of a class of integrable systems in 

OO . Quantum Mechanics. The method is applied in the study of radial oscillator, 

Morse and Coulomb potentials to obtain a wide set of raising and lowering 

^ , operators, and to show clearly the connection that link these systems. 
OO 

m 

O ; I- INTRODUCTION 



Qv . We shall begin this section by recalling some basic facts of the standard factorization 

method, as can be found for instance in |[l|,0], mainly to fix the notation. Afterwards, we 

PL|| will set up the general lines to define more general factorizations, and the way they depart 

from the conventional ones previously characterized. 

j^ ■ Let us consider a sequence of stationary one dimensional Schrodinger equations, labeled 

^ . by an integer number £, written in the form 



H'€ = {-^ + y\^)]€{^) = Ei€{^). (1) 

where the constants % and m have been conveniently reabsorved. If such a set (or 'hierarchy') 
of Hamiltonians can be expressed as 

H' = XtXi - qii) = Xi^.Xl, - g(£-l), (2) 

where 

Xf = T^ + w,{x), (3) 

being w^{x) functions and q{t) constants, then we will say that they admit a factorization. 
From (0) we have that Xf are hermitian conjugated of each other, \X^\ = X/, with 
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respect to the usual inner product of the Schrodinger equation. This is consistent with the 
factorization (|^) and the hermiticity of H^. 

We shall focus our interest in studying the discrete spectrum of each Hamiltonian, so we 
further impose that the equation 

X^4{x)=Q (4) 

will determine the ground state of H^ if it exists. Of course many other properties related 
with the continuous spectrum can also be derived with the help of factorizations, but they 
are out of our present scope. 

Some consequences that can immediately be derived from the previous conditions are 
enumerated below. 

i) Spectrum. Let ip^ be the ground state of H^ as stated in (1), then its energy is 
precisely E^ = —q{i). When there are excited bounded states ijj^, with n = i,i + l, . . ., 
their energy is given by E^ = —q{n). Therefore, in these circunstances, —q{€) should 
be an increasing function on L 

ii) Eigenfunctions. It is straightforward to check, for each ^, the intertwining relations 

H'Xi = X+H'^\ H'^'Xi = XiH'. (5) 

Let us designate by TC^ = {{'ipn}n>i) the Hilbert space spanned by the bounded states 
of H^, for ^ G Z. Then, due to (^ the operators Xf" link these spaces as 

xi : n' -^ n'+' xj : n'^^ -^ W 

(6) 

X^^lix) oc ^^^(x), Xi^l+\x) oc ^i{x). 

Remark that the action of Xf preserve the label ra, that is, they connect eigenfunctions with 
the same energy E"^. If the eigenfunctions are normalized we can be more explicit: up to an 
arbitrary phase factor. 



XtVn^\x) = ^qii)-qin)i;i{x), n>L 

Similar considerations would also apply if the ground states were defined through X^ . 
Depending on each particular problem we will use one of the following notations 

X+i(r)^i, = 0, if ^<0, (8) 

XI,{t)4 = Q, if £>0. (9) 

For such a case —q{tj must be a decreasing function of L We shall also have the opportunity 
to illustrate this situation in some examples along the next sections. 

Now, it is natural to define a set of free-index linear operators {X^, L} acting on the 
direct sum of the Hilbert spaces Ti = ®{H^ by means of 



X-i^l := X^^l X+^l := Xl.^l L^i := £^^, (10) 

where one must have in mind @ and (^. That is, the operators X^ act on each function 
ipn{x) by means of the differential operators (^ changing £ into £ =F 1. The action on any 
other vector of 7i can be obtained from ([T0|) by hnearization, but we shall never need it. 
At this moment we are not in conditions to guarantee that the space 7i is invariant under 
this action (it might happen that the action of X"^ on H could lead us to the continuous 
spectrum, or even to an unphysical eigenfunction), but we postpone this problem to the 
examples of Section III. 

Taking into account our definitions (|TD|), it is straightforward to arrive at the following 



commutators, 

[L, X±] = tX±, [X+, X-] = q{L) - g(L-l). (11) 

It is clear that the set of operators {X^, L} in general does not close a Lie algebra; relations 
(|ll|) only allow us to speak formally of an associative algebra. 

There are many aspects of the conventional factorizations above characterized which can 
be modified, mainly with the objective of being aplicable to a wider class of systems (see for 
example 0). However, in this paper we are interested in going deeply into the possibilities 
of this method on a class of systems where the usual factorization can already be applied, 
so that it could supply us with additional information. With this aim, we shall stress here 
on two points that will be useful in the next sections. 

First, we shall assume that the operators Xf do not have to take necessarily the form 
given in (|^). In particular, if we have a family of invertible operators D^ and define Y^ = 
X/D^^, Yf = DiX^ , we will also have 

H' = XJXi - g(£) = r+r,- - q{£). (12) 

The new factor Di may be a function (which would add nothing specially new) but also a 
local operator, i.e., an operator acting on wavefunctions in the form 

D,i;{x)=i;{g,{x)), (13) 

where gi is a bijective real function. An example of such an operator, which was already 
used in 0, is given by the dilation, 

D{ij,)ijj{x) = ^Pifix), /i > 0. (14) 

Second, an eigenvalue equation can be characterized by more than one label; this con- 
sideration has also been explored by Barut et a/ ||^, but in another context. In the next 
section we shall deal with two real labels; this will enable us to have more possible ways to 
factorize the Hamiltonian hierarchy, and the sequence of labels will not be limited (essen- 
tialy) to the integers, but it will be constituted by a lattice of points in M . This increasing 
of factorizations will reflect itself in a larger algebra of free-index operators. In particular, 
among them, there can be lowering and raising operators for each Hamiltonian, which can 
never be obtained by the conventional factorization method. Section HI will illustrate how 
our general method works when it is applied to three well known potentials: radial oscilla- 
tor, Morse, and radial Coulomb. For each of these potentials we shall see that the results 
so obtained can be used to recover, as special cases, those corresponding to the standard 
factorizations. Finally some comments and remarks will end this paper. 



II. REFINED FACTORIZATIONS 



Once the spectrum E^ of the hierarchy H^ is known, we propose a somewhat more general 
factorization of the eigenvalue equations than that one already displayed in (S), as follows: 



hnA^) 



H'-Et 



BniAn,t - (t){n, i) = A-iB- 



m. 



(15) 



This must be understood as a series of relationships valid for a class of allowed values of the 



parameters (n, 



\ Here Bn/ and An^i are first order differential operators in the wider 



sense specified in the previous section, hn/{x) denote functions, and 4>{n,£) are constants. 
The {n,i) values depend on {n,i), i.e., {n,i) = F{n,i), being F : M^ — >■ M^ an invertible 
map defined on a certain domain. The iterated action of F or F^^ on a fixed initial point 
(no, io) G K originates a sequence of points in M that will play a role similar to the integer 
sequence £ in M for the usual factorizations. In principle the points {n, i) obtained by this 
new approach can take integer values for both arguments, but we do not discard a priori 
other possibilities. 

The problem of finding solutions to this kind of factorizations becomes more involved 
because we have additional functions hn,e{x) to be determined. Nevertheless, an important 
and immediate consequence of (|15D is that the operators Bn,e, An/ share properties similar 
to (^ with respect to their analogs {Xf} : 



KAx){H^ 



Ei 



A. 



n,t 



A 



n,l 



K/{x) (H' - Ei 



5, 



n,l 



where Fin, 



[riy 



h^^x) {h^ - Ei)] = [hnAx) {h'' - E 
. Therefore, using the same notation as in 



B„ 



(16) 



An/ : TC -^ H 

An/i)i{x) (Xlpiix), 



Bn/: n' -^ n' 

Bn/lpiix) CCIpiix). 



(17) 



In this case, the most relevant differences with respect to the usual factorizations are: 



i) Bn/, An/ in general do not preserve the energy eigenvalue, they may change both labels 
n and i. 

ii) An/ does not act on the whole space TC^, it acts just on the eigenfunction ipi{x) G H^ 
(the same can be said of Bn/ with respect to ipHx) G TC^). 

When n = h and hn/{x) = 1, we recover the conventional case with Bn/, An/ playing the 
role of X^ , X'[ , respectively. However, the hermiticity properties for the general case are 
lost because the product Bn/An/ gives not the Hamiltonian operator alone, but it includes 
also a non constant multiplicative factor. 

We can define the free-index operators {A,B,L,N} as we did in (PTJ]), where the latter 
is defined by Nip^ = nipi. They satisfy the following commutation rules 



[L, B] = B{L - L), [N, B] = B{N - N), [B, A] = 0(iV, L) - 0(iV, L) 
[L, A] = {L- L)A, [N, A] = {N- N)A, [N, L] = 0, 



(18) 



where (A^, L) = F{N, L). As the operators L, N commute, their eigenvalues are used to label 
the common eigenfunctions ipHx). We must also notice that the equation An/ipn{x) = (or 
Bn,e'ipi{x) = 0) does not necessarily give an eigenfunction of H^ (or H^); this happens to be 
the case only when 0(n, i) = 0. 

III. APPLICATIONS 

A. Radial Oscillator Potential 

As usual the Hamiltonian of the two dimensional harmonic oscillator includes the effective 
radial potential V^{r) = r^ + ^ 2-2 ~ ; where £ = 0, 1 ... is for the angular momentun. 
The related stationary Schrodinger equation has discrete eigenvalues denoted according to 
the following convention, 



E; = 2n + 2, n = 2u + 



z^ = 0,l, 



It can be factorized in two ways according to our general scheme: 



H'-Ei 



<P , (2f+l)(2£-l) , 



R* 4* — 



n. 



i = 1,2 (19) 



with 0* (n, €) given by 



0i(n,£) = -^(n + £ + 2), 02^n,£) = -l(n-£ + 2), 



(20) 



and where the action on the parameters associated to each factorization is given respectively 
by the functions 



(n, I) = F\n + IJ + 1), (n, £) = F'^{n + 1, £ - 1). 

This can also be written in an easier notation, 

A^ ■ (n, £) -^ (n + 1, £ + 1) A^ ■ (n, £) ^ (n + 1, £ - 1) 

51 : (n + 1, £ + 1) ^ in, £) B^ : {n + l,i - 1) ^ {n, £). 

The explicit form of these intertwining operators is 



(21) 



(22) 



<^(0 = I \i-r-{i+l/2y- 



Bl Ar) 



^+r+(£+l/2)i 



KAr) 



BlAr) 



dr 



-r+(^-l/2)- 



^+,_(£_l/2)i 



(23) 



Observe that in this case, as hn/{r) is a constant, we are able to implement also the her- 
miticity properties {A^^ = —B^. The nonvanishing commutation rules for the free-index 
operators {A^, L, A^, B^; i = 1, 2} are shown to be, in agreement with 



[L, B'] = i-iyB\ [N, B'] = -B' 

[L, A'] = -{-iyA\ [N, A'] = A\ 



[A\ B' 



1,2. 



(24) 



These commutators correspond to two independent boson algebras with A^, L being a hn- 
ear combination of their number operators. Formally we can extend the values of i so to 
include the negative integers. This is physically appealing because in two space dimensions 
(only!) £ represents the L^-component of angular momentum, so that it could take negative 
integer values. Of course the extension ipn^{r) := ipn{r-) above proposed is consistent with 
such an interpretation: (i) The radial components for opposite L^-values have to coincide, 
and (ii) The potential V^ is invariant under the interchange £ — > —L With this conven- 
tion, the Hilbert space Ti of bounded states is invariant under the action of the operators 
{iV, L, A*, 5*; i = 1,2}, so that it constitutes the support for a lowest weight irreducible 



representation for the algebra ( |2^ based on the fudamental state V'^^q- 

It is worth to notice that, taking into account (|2^) , the composition {A^A^^B^B"^} 

constitutes the lowering and raising operators for each Hamiltonian H^ ^ while the pair 

{A^B^, A^B^} connects states of different Hamiltonians H^ with the same energy, changing 

only the label L 

We shall compare briefly the above results with the conventional factorizations of the 

two-dimensional radial oscillator potential 0. It is well known that there are two such 

factorizations which we will write in the form: 

(a) X+X,~ - g,(£) = Hi = H' - 2i (25) 

(6) Z+Zi-q,{i)=Hi = H' + 2i, (26) 

with H^ = — ^+r^-|- ^ "'"2.2 ~ • Then we have the following identification: 
Case (a) 

1. Operators: X/ = -2Bl^, X^ = 2v4^_^, q^{i) =Ai-2. 

2. Ground states: X^"^iV^i^ = 0, £ < 0. 

3. Energy eigenvalues: E^ = An + 2, with n G Z^ and n > —L 
In this case we have used a notation in agreement with (H). 
Case (6) 

1. Operators: Z+ = -2Al_^^^^^, Zj = 2Bl_^^^^^, q,{i) = -U~2. 

2. Ground states: Z^ipl = 0, £ > 0. 

3. Energy eigenvalues: El^ = An + 2, with n G Z^ and n> L 

Therefore, as there is a correspondence between the results of the conventional and our 
factorizations, one might conclude the total equivalence of both treatments. However, we 
make a remark worth to take into account: the conventional factorizations make use of 
two Hamiltonian hierarchies, H^ and iff, whose terms differ in a constant 4£, while the new 
factorizations use only one H^. If we want that both factorizations (a) and ih) be valid inside 
the same hierarchy it is necessary to adopt the properties of our approach in the following 
sense: either the operators Xf or Zf (or both pairs) must change not only the quantum 
number I but also n. In this way we have shown, by means of this simple example, that 
the factorizations presented here prove to be quite useful providing directly a more natural 
viewpoint. 



B. Morse Potential 

In this case we have eigenvalue Schrodinger equations for the whole real line x G M with 
the potentials 

V\x) = (^) (e^"^ - 2{i + 1) e"^) , a>0,i>0. (27) 

Often in the literature the Morse potentials are written V{y) = A (e^^"^ — 2e^"^). This 
form can be reached from ( [27| ) by a simple change of the variable x = —y+k, with e"'^ = £+1. 
The energy eigenvalues can be expressed as 

2 

E^ = -^n^ n = £-2u>0; z/ = 0,l,2... (28) 

In order to have bounded states it is necessary the restriction £ > 0; the critical value £ = 
has in this respect an special limiting character, and it is convenient to take it into account 



as we shall see later. According to (pS]), the eigenfunctions ip^ are characterized by labels 



satisfying n < £; this means that the ground states will be defined through 
There are two new factorizations 



a^ 



^ -K 



Bl/x)Al,ix)-<P\n,£), 2 = 1,2, (29) 



with (f)^{n,£) given by 



<P\n,£) = -^{£ + n + 2), ^\n, £) = -^{£ ~ n + 2) , (30) 

and the action on the parameters (n, £) for each factorization by the functions 

{n,£) = F\n + l,£ + l), {n,£) = F\n - 1,£ + 1). (31) 

The explicit form of the intertwining operators (^) is 

Rl (t^ — 6-""=/^ d I 1 ax/2 I n+1 -ax/2 
^n/\-^) — a dx^2^ ^ 2 ^ 



Al (^:\ _ e-°^/^ d l^ax/2 n^-ax/2 
^nAX) - — ^d^-2^ ~2^ 

d2 (^\ _ e-"^/^ d ,l^ax/2 n-l ^-ax/2 
^nA^) - — ^d^+2^ 2~^ 

42 ( \ _ e-°^/^ d _1 ax/2, n -ax/2 
^nA^) - —^d^ 2^ +2^ 



(32) 



(33) 



As in the oscillator case we have two pairs of operators that change simultaneously two 
types of labels: one, £, is related to the intensity of the potential, although here it can not 
be interpreted as due to a centrifugal term. The second one, n, is directly related to the 
energy through formula ( ^8]) . The (nonvanishing) commutators of the free-index operators 
are 

[L, B^] = ~B\ [N, B'] = i-iyB\ [A\ B'] = 1 

[L,A^] = A\ [N,A'] = ~{-iyA\ i = l,2. 



Observe that in this case the function hn/{x) = — e~"^/a^ is not a constant, so the 
hermiticity relations among the operators {A^, B^; i = 1,2} are spoiled. Let us take i G Z"*", 
and formally allow for negative n-values in (^), i.e., ±n = i — 2i>; this is admissible because 
in the operators of (^)-(|33|) we have a symmetry under the change n — >■ —n. Then the 



Hilbert space TC of bounded states enlarged with the (not square integrable) states V'n=oi 
i = 0,1,2 .. ., will be invariant under the action of all the operators defined in this section. 
The lowest weight state is played in this case by a not square-integrable wavefunction, ip^Q. 
We can of course build other operators out of the previous ones, changing exclusively 
one of the labels: the pair {A^A'^,B^B^} change i (in +2 or —2 units, respectively), while 
{A^B'^,A^B^} change n (also in +2 or —2 units, respectively). It is interesting to show 
explicitly the form taken by the former couple: 

where {A^A^)n,e = Al^_^^^^Af^^ and {B^B'^)n/ = -B^^-B^^.^^.,.]^, according to the rules of the 
action of free index operators (|T3), (pil|). They can be identified with the usual factorization 
operators for the Morse Hamiltonians H^ described in the first section in the following way 

1. Factorization: X+X^ - q{i') = H^^' , f G Z+. 

2. Operators: X+ = -a iB'BX,2i' , X,7 = -a {A'AX,2i' , q{i') = «'(f + 1)' • 

3. Ground states: X^_^tpf, = 0, f > 0. 

4. Energy eigenvalues: E2^i = —a'^{n'y, with n = 2n', n' G Z^, and < n' < C' . 
This time the notation, as it was mentioned above, is in agreement with (|^). 

C. Radial Coulomb Potential 

After the separation of the angular variables, the stationary radial Schrodinger equation 
for the Coulomb potential in two dimensions takes the form 

H'€{r) = 1-^ + ^'^^yf"'^ - - W^(r) = K€{rl (36) 



dr"^ Ar"^ 



r 



where the values of the orbital angular momentum are positive integers £ = 0,1,2... 

The computation of the discrete spectrum associated to the bounded states of H^ can 
be easily obtained by means of the conventional factorizations (^ with 

..+ d 2£ + l 2 ,,, -1 

X^ = T-, + —. , qU) = 9 • 37) 

^ ^dr 2r 2i + V ^^ ' {i+l/2f ^ ' 

Therefore, according to the results quoted in Section I, we have 



When our method is apphed to the hydrogen Hamiltonians H^ of equation (|36D with the 



eigenvalues E^^ (|38|) , we obtain two independent solutions that read as follows 

(2n + l)r 



BitKi + ^ + n + l 



H'-Et 



r2 42 



+ n 



(2n + l)r r 



H'~Ei 



(39) 
(40) 



The explicit form of the operators {A\ i?*}j=i,2, is displayed below: 

Bi,{r) = (2n + 1)^^ (4^|:+^+^) c^-V^ Z)(c„) 



5^,(r) = (2n + 1)^^ (4^i+£^-^) 0'^^ D{c„ 



Aleir) = Dic^-^)ci/H2n + lY/'{r^i 



___rVi I 2^+1 
dr 2n+l~'~4rl/2 



The symbol -D(/i) in (0)-(H3) is for the dilation operator (|14|), and c„ 



2n+2 
2n+l' 



(41) 



(42) 



Thus, in 



this example we are dealing with general first order differential operators as explained in 
Section I. For the first couple {A^,B^} we have (n, £) = (n + 1/2, £ + 1/2), while for the 
second pair {A^ , B^}, (n, i) = {n + 1/2, i - 1/2). 

The nonvanishing commutators among the free-index operators are 



[N,B' 
[N, A\ 



2 

2/1 , 



[L, B^] 
[L, A^] 



{-iylB\ 

-i-iy^A^ 



[A\B^]=I, 
i = 1,2 . 



(43) 



In other words, as in the previous examples, we have a set of two independent boson oper- 
ator algebras. The problem with these operators is that they change the quantum numbers 
(n, i) in half-units, so that they do not keep inside the sector of physical wavefunctions. To 
avoid this problem we can build quadratic operators {A^A^,B^A^,B^B^ij=i2 satisfy- 
ing this requirement; such second-order operators close the Lie algebra sp(4, M) 0, which 
includes the subalgebra su{2) (whose generators connect eigenstates with the same energy 
but different £'s). It is worth to write down these quadratic operators: 



{B'A%j^ 
{B'AX,i 



(2^+l)(2n+l) (id. 2^+1 ]_ 

2 \2dr~^ 4r 2^+1 



(2£+l)(2n+l) 
2 



V 2dr^ 



2i+l 



Ar 



2£+l 



(44) 



They constitute, up to global constants, the usual factorization operators given in 
X/ oc {A^B'^)n/, X^ oc {A'^B^)n/. Another subalgebra is su{l, 1) (relating states with the 
same i but different energies or n values). Once included the negative i values, as we did for 
the radial oscillator potential, the space Ti is the support for what it is called a 'singleton 
representation' [^ of so(3, 2) ^ sp(4, R). There is one lowest weight eigenvector V'n=o ^ ^^ 
from which all the representation space is generated by applying raising operators. 



IV. CONCLUSIONS AND REMARKS 

We have shown that a refinement of the factorization method allows us to study the max- 
imum of relations among the Hamiltonian hierarchies that the conventional factorizations 
are not able to appreciate. The operators involved obbey commutation rules which show 
the connection existing among the three examples dealt with in this paper: they have the 
same underlying Lie algebra associated with confiuent hypergeometric functions. In other 
occasions the conventional factorizations have been used in this respect, but we have seen 
that such an approach is partial and not complete at all. 

Usually the Hamiltonian hierarchies are obtained from higher dimensional systems after 
separation of variables (or by any other way of reduction). Such systems have symmetries 
that are responsible for their analytical treatment. These symmetries are refiected in the 
many factorizations that the hierarchies can give rise to by means of the thecnique we have 
developed. We have limited our study to A^ = 2 space dimensions for the radial oscillator and 
Coulomb potentials because they are the simplest cases to deal with. For other dimensions 
there appear certain subtleties, in the sense that the Hilbert space 7i of bounded states is 



no longer invariant under the involved operators [|Il| . 

Finally, let us mention that we have limited ourselves to some examples (all of them 
inside the class of shape invariant potentials |T^), but it is clear that the whole treatment 
is aplicable to the remaining Hamiltonians in the classification of Infeld and Hull p[ . 
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